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We propose a general procedure for implementing dynamical decoupling without requiring arbitrarily strong, 
impulsive control actions. This is accomplished by designing continuous decoupling propagators according to 
Eulerian paths in the decoupling group for the system. Such Eulerian decoupling schemes offer two important 
advantages over their impulsive counterparts: they are able to enforce the same dynamical symmetrization but 
with more realistic control resources and, at the same time, they are intrinsically tolerant against a large class of 
systematic implementation errors. 
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PACS numbers: 03.67.-a, 02.70.-c, 03.65.Yz, 89.70,+c 

Dynamical decoupling provides a well-defined framework 
for addressing a variety of issues associated with the manipu- 
lation of open quantum systems and interacting quantum sub- 
systems. Inspired by coherent averaging methods in nuclear 
magnetic resonance spectroscopy [fij], and cast in control- 
theoretic terms in [|[ |], decoupling techniques are attract- 
ing growing interest from the quantum control and quantum 
information processing (QIP) communities. Significant appli- 
cations have resulted in the area of reliable QIP, where de- 
coupling has been instrumental in the development of quan- 
tum error suppression schemes [|| B], with the potential for 
noise-tolerant universal quantum computation on dynamically 
generated noiseless subsystems |Q. In addition, variants of 
the basic decoupling concepts play a role in protocols for uni- 
versal quantum simulation of both closed- and open-system 
dynamic s M ||, ffl, with implications for encoded simula- 
tion [[ToJ |1 In a broader context, applications of dynamical 
decoupling to problems that range from inhibiting the decay 
of unstable states [O], to suppressing magnetic state decoher- 
ence [n3J, or reducing heating effects in linear ion traps [|l4j] 
have been recently envisaged. 

From the point of view of implementation, dynamical de- 
coupling has relied on the ability of effecting sequences of 
arbitrarily strong, instantaneous control pulses. That is, it re- 
quired the ability to impulsively apply a set of control Hamil- 
tonians with unbounded strength (the bang-bang (b.b.) as- 
sumption While providing a convenient starting point, 
such a scenario suffers from being extremely unrealistic for 
applications. In a physical control setting, additional disad- 
vantages associated with b.b. decoupling include the diffi- 
culty of simultaneously describing the evolution under the 
natural (drift) Hamiltonian and the control terms, as well as 
the poor spectral selectivity of b.b. pulses, with substantial 
off-resonance effects. Finally, although compensation tech- 
niques based on composite rotations exist for stabilizing con- 
trol pulses against operational imperfections Jl5|], they are 
hard to reconcile with the b.b. framework, which does not 
easily lend itself to incorporating robustness features. 

In this Letter, we overcome the shortcomings of the b.b. 
formulation by showing how to implement dynamical decou- 
pling based on continuous modulation of bounded-strength 
Hamiltonians. If Q is the discrete group specifying the desired 



b.b. decoupler, the basic idea is to constrain the motion of the 
control propagator during each cycle along a path that inter- 
polates between the elements of Q. Under mild assumptions 
on the control Hamiltonians, a decoupling prescription induc- 
ing the same symmetry structure as in the b.b. limit can be 
constructed by exploiting Eulerian cycles on a Cayley graph 
of Q. In addition to significantly weakening the relevant im- 
plementation requirements, Eulerian decoupling turns out to 
be largely insensitive to control faults, opening the way to the 
robust dynamical generation of noise-protected subsystems. 

Decoupling setting.— Let the target system S be defined 
on a finite-dimensional state space Hs, and let End(Hs) 
be the corresponding operator algebra. Thus, Hs — C d , 
End(Hs) — Matd(C) for some d, with d = 2™ for an n-qubit 
system. S may be coupled to an uncontrollable environment 
E, whereby the evolution on the joint state space Hs ®He 
is ruled by a total drift Hamiltonian Hq = Hs (E> He + Is <8> 
He + J2 a Sa ® E a for appropriate traceless noise genera- 
tors S a 6 End(Hs) ®- A decoupling problem is concerned 
with characterizing the effective evolutions that can be gener- 
ated from Ho via the application of a control field H c (t) Cg> 
acting on S alone [pp. Let the control propagator be 



U c (t) = Texpj-i^ dt' H c (t')} 



(1) 



with h = 1. In a frame that removes the control field, the dy- 
namics is governed by a time-dependent Hamiltonian H(t)= 
Ul(t)HoU c (i), and the overall evolution in the Schrodinger 
picture results from the net propagator 



U(t) = U c (t)Texp{-i f dt'H(t')} 



(2) 



Assuming that the control action is cyclic, U c (t + T c ) = U c (t) 
for some cycle time T c > and for all t, the stroboscopic 
dynamics U(t M = MT C ), M £ N, can be identified with 
the effective evolution induced by H(t) in (^J). First-order 
decoupling aims at generating the desired evolution to lowest 

order in T c , U(tM) = cxp(—iH tm), where 



H ia) = -\ dt'Ul(t')H U c {t') 



(3) 



2 



While higher-order corrections can be systematically evalu- 
ated, the approximation (|3|) tends to become exact as the fast 
control limit T c — > is approached |[ |[ ■ 

In the simplest b.b. decoupling setting, the time-average in 
(H) maps directly into a group-theoretical average. Let Q be a 
discrete group of order \Q\ > 1, Q = {gj}, j = 0, . . . , \Q\ — 1, 
acting on Hs via a faithful, unitary, projective representation 
/i, C U{TLs)- Let images of abstract quantities under fi 
be denoted as fJ,(gj) = gj, and so forth JT^]. Then b.b. decou- 
pling according to C/ is implemented by specifying U c (t) over 
each of the \Q\ equal sub-intervals defining a control cycle [fj|j: 



f 7 , 



l)At 



€[0,At), (4) 



with T c = |£| At for At > 0, and I = 1, . . . , The result- 
ing control action corresponds to extracting the (/-invariant 

component of Hq, — YL g (Ho), where 



U G W = FT E a] X g 3 , X G End(H s ) , 



(5) 



is the projector onto the commutant CQ' of C£ in End(Hs) [j^, 
0], According to (Q), J7 c (t) jumps from to gi = 

{geg\-i)ge-i through the application of an arbitrarily strong, 
instantaneous control kick at the ^'th endpoint tp = tAt, real- 
izing the b.b. pulse pp = geg\_i [§]■ 

Eulerian dynamical decoupling.— We seek a way for 
smoothly steering U c (t) from gp-\ to gp by a control action 
distributed along the whole ^'th sub-interval. Let T = {7a}, 
A = l,...,|r|bea generating set for Q. The Cay ley graph 
G(Q,T) of Q with respect to Y is the directed multigraph 
whose edges are coloured with the generators jl7|, where ver- 
tex gi-\ is joined to vertex gi by an edge of colour A if and 
only if gpgj\ = 7a i.e., g e = i\gt-i. Physically, imagine 
that we have the ability to implement each generator j\, by 
the application of control Hamiltonians h\(t) over At, 



7A 



Tcxp{-i J dt'/i A (t')} , 



A = 1, 



(6) 



The choice of h\{t) is not unique, allowing for additional 
implementation flexibility. Once a choice is made, the con- 
trol action is determined by assigning a cycle time and a rule 
for switching the Hamiltonians h\(t) during the cycle sub- 
intervals. We show how a useful rule results from sequen- 
tially implementing generators so that they follow a Eulerian 
cycle on G(Q, Y). A Eulerian cycle is defined as a cycle that 
uses each edge exactly once $1%. Because a Cayley graph is 
regular, it always possesses Eulerian cycles, whose length is 
necessarily L = \G\\Y\ [17|. 



Let a Eulerian cycle beginning at the identity go of G be 
given by the sequence of edge colours used, Ve = {pe)e, 
with £ = 1, . . . , L, and pp = 7a for some A, for every t. 
Note that each vertex has exactly one departing edge of each 
colour, so that Ve determines a well defined path. We define 



Eulerian decoupling according to Q by letting T c = LAt and 
by assigning U c (t) as follows: 



U, 



l)Af + s 



ut{s)U c ({t-l)At) , 



(7) 



where s e [0, At), and ue(s) = Tcxp(-i J Q S dt' hp(t')), 
up(At) = pe, I = 1, . . . , L. This decoupling prescription 
means that during the £'th sub-interval one chooses as a con- 
trol Hamiltonian the one that implements the generator j\, 
with 7a colouring the edge pp in Ve- The effective Hamilto- 
nian H ^ under Eulerian decoupling is obtained by evaluating 
the time average (fy based on (Q). The resulting L terms can 
be partitioned into |T| families, each corresponding to a fixed 
generator 7a. Because Ve contains exactly one 7A-coloured 
edge ending at any given vertex gj , each family effects a sum 
over the group elements as in Thus, the quantum opera- 
tion Qg corresponding to (Q) can be expressed as 



Q d (X)=IL d (F f (X)), XeEnd(H s ), 



(8) 



with the map Fp implementing an average over both the group 
generators and control sub-interval: 



ir 



i r 1 i r At 

F r( X ) = TF\Y, A - t dsu{(s)Xu x (s). (9) 

' ' A=l J ° 



Thanks to the way enters (§), Q§(X)e CQ' for an 
arbitrary input X. This property will be repeatedly used in 
the following. The link between Eulerian decoupling and Q- 
symmetrization is established upon enforcing some additional 
compatibility between and Fp. 

Theorem. Let X be any (time-independent) operator on 
7"ts, and let be defined as above. If the controls are 

chosen in the decoupling group algebra, hp(t) G <CQ for all 
t e [0, At] and for all £ = 1,...L, then 

Q g {X) = Yl g (X) , X e End(Hs) . 

Proof. The assumption on the controls implies that U\(s) E 
CG VA, Vs e [0, At]. Thus, F f (Y) = Y for every (time- 
independent) Y G CQ'. Now let X G End(?is) an( J calculate 

Q 2 § (X) = Yl g (Q g (X)) = Yl g (F f (X)) = Q g (X). Thus, Q g 

is a projector. Because Range Qq C CQ' , Q g = YLz iff Q g 

has identity action on CQ'. Let Y G CQ', then Q g (Y) = 
Yl g (F f (Y))=Yl g (Y). M 
The b.b. limit is formally recovered by letting Ff be the 
identity map. In the Eulerian approach, at the expense of 
lengthening the control cycle by a factor of the same Q- 
symmetrization can be attained using bounded controls. The 
maximum strengths achievable in implementing the genera- 
tors (||) directly affects the minimum attainable T c , and there- 
fore the accuracy of the averaging [Bp. While the overhead |T| 
depends on the specific group, it is worth noting that, similar 
to {J, Q g satisfies the property that Qq(X) = Qgj^^ (X) 
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whenever Qq is a normal subgroup of Q and A e CQo' JT^]. 
Thus, if the dynamics is already C?o -invariant, Eulerian decou- 
pling according to Q can be accomplished by using a Cayley 
graph of the smaller quotient group G/Go ■ 

Robustness analysis.- The fact that control actions are now 
distributed along finite time intervals translates into major 
gains in terms of resilience of Eulerian schemes against imper- 
fections in the controls themselves. Imagine that systematic 
implementation errors result in a faulty control Hamiltonian 
H' c (t), and partition H' c (t) into 

H' c (t) = H c (t) + AH c (t) , (10) 

such that H c (t) (zCQ is the intended control Hamiltonian, and 
AH c (t) is the error component. Now work in the same frame 
used earlier, which only removes the ideal control part from 
the effective Hamiltonian. Because H(t) = Hq + H' c (t) = 
[Ho+ AH c (t)] +H c {t), this maps the evolution under Hq with 
the faulty control H' c (t) into the evolution under Ho + AH c (t) 
with the ideal control. Thus, the new effective dynamics is 
obtained by replacing Hq with Hq + AH c (t) in (^|). 

Suppose that the faults are properly correlated with the un- 
derlying path, meaning that every time a particular generator 
7a is implemented, the same imperfection occurs at equiva- 
lent temporal locations within the sub-interval, regardless of 
the position of 7> along Ve- Then AH C ((£ — l)At + s) = 
Ah\(s), A being the colour of the edge that Ve uses during 
the £'th sub-interval. By a similar calculation as in the ideal 
case, the quantum operation Qg is modified as follows: 

Q' g (X) = U g (X) + Q Q {AH C ) , AeEnd(H s ), (11) 

where Qg(AH c ) can be computed as in (||) and (^|), but with 
the operator X in the integral replaced by one that depends 
on s and A. Thus, Qg(AH c ) is a functional of the fault his- 
tory over [0, At], which characterizes the residual control er- 
rors experienced by the system. Notably, two useful features 
emerge: without extra assumptions, such residual control er- 
rors belong to CQ'. If, in addition, AH c (t) is itself (as H c (t)) 
in CQ, then all control effects remain in CQ, and the residual 
control errors belong to the center Z^~„ = CQ n CQ'. 

The effects of Qg(AH c ) may still adversely impact the per- 
formance of the system. However, they can be compensated 
for by encodings in appropriate subsystems [g]. Let J g J 
label the irreducible components of CQ. Then Tts can be rep- 
resented as 

Us (BjHj ~ 0j Cj ®Vj ~ ©j C nj ® C dj , (12) 

with nj, dj e n jdj = d, and the action of the de- 

coupling group algebra and its commutant given by CQ ~ 
®jt nj ® Mat d/ (C), CQ' ~ ©jMat n/ (C) <Z> t dj , respec- 
tively. Because both Hg(S a ) and Qg(AH c ) are in CQ', T>j- 
subsystems are noiseless and their dynamical generation ro- 
bust regardless of whether AH c (t) belongs to CQ or not. This 
applies in particular if Q acts irreducibly on Hs, in which case 



a robust implementation of maximal decoupling is achievable 
by averaging over a nice error basis on C d [||, 0. In fact, en- 
coding into 2?j-subsystems may be valuable even in situations 
where the assumption that the controls are in CQ cannot be 
met: as Qg(S a ) £ CQ', T> j-subsystems remain unaffected by 
the noise. Note that for such subsystems, both the implemen- 
tation of the decoupling scheme and the execution of encoded 
control operations are to be effected through fast modulation 
of Hamiltonians along the control cycle ||[ [§. 

Whenever Q^(AH C ) originates from faults in CQ, addi- 
tional options are viable. If the representation fi is primary, 
Z^, = CI, then any systematic error is effectively elim- 
inated, and no encoding is necessary as long as noise sup- 
pression is ensured, that is, Hg*(S a ) = for all a. If /x is 
not primary, then elements in the center are diagonal over 
each irreducible component. Thus, encodings into either a 
Ti. j-subspace or a Cj-subsystem are insensitive to the control 
faults and protected against the noise generator if Hg(S a ) € 
Z^-g as well. In practice, choosing a Cj-subsystem may be es- 
pecially appealing, because not only is universal encoded con- 
trol achievable by less-demanding, slow application of Hamil- 
tonians in C^^0|, but added robustness against arbitrary con- 
trol errors in CQ is automatically provided [^]. Next, we out- 
line some applications relevant to QIP. 

Example 1: Eulerian Carr-Purcell decoupling on a qubit.- 
Consider a single decohering qubit, {S a } = {u z } [§]• The 
decoupling group Q = = {0, 1} is represented in U(C 2 ) 
as Q = {l,a x }. There is one generator, 71 = 1, hence 
L = 2 with no overhead with respect to the b.b. case. Let 
u x (s) =Texp(— i J Q dt' h x (t')), for a Hamiltonian h x (t) £ 

CQ realizing 71 = a x (§). If h x (t) = f{t)a x for some f(t), 
any choice such that | J dsf(s)\= ir/2 is acceptable. On 
G(Q, T) choose Ve = (71, 7i)- Then Eulerian decoupling is 
accomplished by letting U c (t) = u x (t), for t £ [0, At), and 
U c {t) = u x {s)a x fort 6 [At, At + s), s g [0, At). By 
explicit calculation of Qg(AH c ), one sees that systematic er- 
rors along cr y , cr z produce no effect. Elimination of residual 
control errors in Z^~ Q requires using the full Pauli group. 

Example 2: Eulerian Pauli decoupling on qubits.- Let 
Q = {1, X, Y, Z} be the Pauli error basis for a qubit, with 
X = a x , Z = a z , and Y = XZ. This corresponds to 
Q = Z2XZ2, projectively represented in U(C 2 ). Q has 
two generators, e.g. 71 = (0,1), 72 = (1,0), realized as 
71 = A", 72 = Z, respectively. An Eulerian path on G(G, T) 
isV E = (71,72,71,72,72,71,72,71)' of length L = 8. The 
assumption that both h\ and Ah\, A = 1, 2, are in CQ is auto- 
matically satisfied, as CQ = Mat2(C). Then (Q) results into a 
robust implementation of maximal averaging, Hg(a u ) = 0, 
u = x,y,z. For n qubits, Q = Z^xZ^, with d = 2™. 
Thus, \Q\ = 4™ and since 2 generators are needed for each 
qubit, L = n2 2n+1 , causing the procedure to be (as in the b.b. 
limit [|]|) inefficient. 

Example 3: Eulerian collective spin-flip decoupling.- For 
n qubits, let Q = Z2X.Z2 act via the n-fold tensor power 
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representation in U((C 2 ) 
and regular for n even. 



where X 



which is projective for n odd, 
For any n, Q = {1, X, Y, Z}, 



and Y = XZ. De- 



coupling according to Q averages out arbitrary linear noise, 



n 



g(S a ) — 0, S a G span{(7„ } |Q. For Eulerian implemen- 
tation, the same path of Example 2 may be used, under the 
appropriate realization of the collective generators 71 = X, 
72 = Z. Ensuring (y-symmetrization requires that the control 
Hamiltonians hi^{t) G CQ. Because both CQ' and Z^g are 
non-trivial, residual control errors arise due to Qg (AH C ). The 

situation is simpler for n even, as CQ is abelian hence support- 
ing four (71 — 2)-dimensional irreducible subspaces Hj. Be- 
sides being noiseless in the decoupling limit and insensitive to 
arbitrary control errors in CQ, encoding into a Hj-subspace 
is further motivated by the possibility to achieve encoded uni- 
versality via slow application of two-body Hamiltonians in 
CQ' [01. For n odd, both Cj and Vj factors may occur. Leav- 
ing aside details here, we note that Dj -subsystems may be 
useful if implementing 71 , 72 via Hamiltonians in CQ is diffi- 
cult in practice. 

Example 4: Eulerian symmetric decoupling.- Let Q = S n 
be the symmetric group of order n, acting on Hs — (C 2 )®" 

via 9j®k=i iV'fe) = ®k=Mg j (k))' 9j S S n . In particu- 
lar, the action corresponding to a transposition (k — 1 k), 
k G {l,...,ft}, effects an exchange gate between qubits 
k — l,k, denoted by swapfc_i.fc. Symmetric decoupling al- 
lows, in principle, to engineer collective error models on S 
starting from arbitrary linear interactions between S and E 
[Q, 0J. A minimal generating set for S„ is given by 71 = (1 2), 
72 = (1 2 ... n) i.e., an adjacent transposition and the cyclic 
shift, respectively. CS n contains the Heisenberg couplings 
h(k,l) — <7fc • 07. In fact, every operator in CS n can be 
realized by applying Heisenberg Hamiltonians [p^|. Focus, 
for instance, on iS3-symmetrization, which may be relevant 
for inducing collective decoherence on blocks of 3 qubits 
[^0|]. Then 71 = swapi^ and 72 =swapi j 2 swap2.3, with 
L = 12. Because exp(— iirh(k, I) /A) =swap / t i z, 71 can 
be implemented by choosing h\ = a\h{\,2), with strength 
ai = 7r/4A<, while 72 can be realized by a piecewise- 
constant Hamiltonian h 2 (t) = a 2 h(2, 3) for t G [0, At/2), 
h 2 (t) = a 2 h(l, 2) for t G [At/2, At], a 2 = 7r/2At A Eule- 
rian path on G(S 3 ,T) is V E = (72,72, 72,71, 72,71, 7i,72, 
7i, 7i, 72, 7i )■ Eulerian decoupling (Q) then allows for a ro- 
bust dynamical generation of the smallest non-trivial noiseless 
subsystem [ pTj p2[ |, supported by a factor Vj ~ C 2 carrying 
the two-dimensional irreducible component J = [2 1] of 53. 

Conclusion.- We developed an approach to dynamical de- 
coupling that combines the group-theoretical essence of the 



b.b. setting with graph-theoretical control design according to 
Eulerian cycles. Besides allowing for considerable leeway in 
the physical implementation of the basic control generators, 
Eulerian decoupling eliminates the need for unfeasible b.b. 
pulses and naturally incorporates robustness against realistic 
control faults. Combined with quantum coding techniques, 
our results significantly improve the prospects that dynamical 
decoupling becomes a practical tool for reliably controlling 
quantum systems and quantum information. 
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